Directing Brownian motion on a periodic surface 
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We consider an overdamped Brownian particle, exposed to a two-dimensional, square lattice 
potential and a rectangular ac-drive. Depending on the driving amplitude, the linear response to 
a weak dc-force along a lattice symmetry axis consist in a mobility in basically any direction. In 
particular, motion exactly opposite to the applied dc-force may arise. Upon changing the angle of 
the dc-force relatively to the square lattice, the particle motion remains predominantly opposite to 
the dc-force. The basic physical mechanism consists in a spontaneous symmetry breaking of the 
unbiased deterministic particle dynamics. 

PACS numbers: 05.45.-a, 05.60.-k, 05.40.-a 



Brownian particle dynamics in two-dimensional peri- 
odic potential landscapes arise in a large variety of dif- 
ferent contexts. Examples include driven vortex lattices 
[11, S S] , surface diffusion [3| , a ring of several Josephson 
junctions colloidal particles or globular DNA in struc- 
tured microfluidic devices [1, 0| and in optical [1, Q or 
magnetic [loj lattices, enzymatic reaction cycles driving 
molecular motors [Tl| . nanoscale friction [T3 | and super- 
lubricity They have recently attracted considerable 
theoretical [l|, and experimental [1| interest for par- 
ticle sorting in two-dimensional periodic structures with 
the help of an externally applied dc-force, whose angle 
relatively to the periodic potential can be parametrically 
changed. The key point is that the resulting particle ve- 
locity may exhibit a different direction than the applied 
dc-force and that the deflection angle may be different 
for different particle species. While the deflection angles 
between force and velocity remain bounded to relatively 
small values, our present system will lead to (practically) 
arbitrary deflection angles. 

A second recent series of papers 0, considers the 
same system but in the presence of an additional circular 
ac-drive. Deflection angles up to 90° (absolute trans- 
verse mobility) have been found in 0], while [lo| reports 
transporting orbits in the absence of a dc-drive. A related 
variant is to replace the circular by a more common, lin- 
ear ac-drive, but now breaking the time-space symmetry 
by chosing a bi- harmonic driving signal and focus- 
ing on the low friction regime [l5|. The setup we will 
consider here is related but simpler: a standard ac-drive 
without any concomitant space-time symmetry breaking 
and negligible inertia effects. 

Unbiased far from equilibrium dynamics of a Brownian 
particle, responding to a dc-force by a directed trans- 
port opposite to that force, have been extensively in- 
vestigated under the label "absolute negative mobility" 
(ANM) 0, Our present work represents the natu- 
ral extension, namely an unbiased far from equilibrium 
system admitting an easily controllable linear response 
into (practically) any direction relative to the dc-force, 
including ANM as a special case. 

We consider the following 2d Langevin dynamics of a 




FIG. 1: Schematic illustration of our model ([T}. The contour 
lines represent the potential U (r) with a cut-off for better 
visualization (white discs) . The double arrow indicates the ac- 
drive A{t)ea, one arrow the dc-bias FStf,, and another arrow 
the particle velocity ([2]). The particle is sketched by the black 
dot. The angle ij) quantifies the "deflection" of v from Fe^. 



Brownian particle with coordinates r 



fit) ^ A(t)e„ + Fe^ - VC/(f(t)) + V2r^{t) 



(1) 



Thus, inertia effects are neglected (overdamped dynam- 
ics) , and the friction coefficient is absorbed into the time 
unit. As illustrated with Fig. [H A{t) is the ac-driving 
signal along the direction := (cos a, sin a), and anal- 
ogously for the dc-bias Fe^. The periodic potential is 
represented by U{r} and thermal fluctuations of temper- 
ature r are modeled by the two delta-correlated, Gaus- 
sian noise components of ^{t). We focus on the particu- 
larly simple rectangular driving A{t) ~ a sign{cos(il i)} 
with amplitude a and period T ~ 2t:/^1. We veri- 
fied that a sinusoidal A{t) leaves all our main findings 
qualitatively unchanged and expect the same for even 
more general A{t). Further, we focus on the potential 
^(^) ~ '^U{f+ Lne-x + LmCy) with r j?^, U{f) — 
u exp{— r/A}/r, and m, A > 0. In other words, we con- 
sider a square lattice of repulsive Yukawa potentials, a 
standard choice for screened charges Again, we ex- 
pect that our results remain qualitatively unchanged for 



-0.2 0.2 0.4 0.2 0.4 

V, [L/T] V, [L/T] 

FIG. 2; Velocity ([2} from numerical simulations of ([1} for 
a = 0.15 ■ 360° = 54°, a G [4.3, 8], = 0°, F = 0.03, n = 4.3, 
r = 2.2 • 10"" (left) and T = 6.7 • 10"'' (right). Shown is 
the parametric dependence of the velocity components Vx, Vy 
on the ac-amplitude a in units of L/T {L = 1, T = 2n/Q). 
Arrows indicate increasing a-values. Upon further decreasing 
-F, a close to linear response behavior of v results (not shown). 

more general U{r), modelling e.g. pinned vortices [l|, 
optical tweezers [1, |1] , or magnetic bubbles [l^ , and we 
have explicitly verified this for Gaussian shaped repul- 
sive and attractive potentials. Without loss of generality 
we choose length and time units with L = 1 and u = 1. 
Regarding A, we obtained practically indistinguishable 
results for all A > 4, variations by a few percent down to 
A ~ 1, and notable quantitative but no qualitative differ- 
ences at least down to A ~ 0.1. In the following wc focus 
on the representative example A = 4. 

The observable of main interest will be the time- 
averaged particle velocity 

I r°° ^ 

V = Vxe-x + VyGy :~ lim - / dt' r{t') . (2) 

t^OO t Jq 

being independent of the seed r(0) and the realization of 
in © for any F > due to ergodicity reasons. 

Generally speaking, the periodic potential, the ac-drivc 
and the dc-bias in ([T]) give rise to several "competing di- 
rections" , whose net effect on the velocity v from ^ is far 
from obvious. For zero bias F, the ac-forcing still keeps 
the system off equilibrium but any non-zero velocity v is 
prohibited by symmetries [l7| . Our first objective is the 
linear response of to a weak dc-bias along the x-axis, cf. 
Fig. [T] Our findings in Fig. [2] exhibit a quite intriguing 
behavior. Keeping all "competing directions" fixed and 
solely changing the ac-amplitude by 40%, almost any di- 
rection of V may arise, even motion exactly opposite to 
the apphcd dc-bias (ANM). 

Next, we focus on a set of parameters close to the oc- 
currence of ANM in Fig. [5] and now ask for the response 
of V upon changing the direction of the dc-bias. Again, 
the results in Fig. [3] arc quite non-trivial, the most re- 
markable feature being that the projection of the veloc- 
ity along the dc-bias is mostly negative, i.e. the particle 
motion remains predominantly opposite to the dc-force 
(90° < V < 270°). Similarly as in Fig. d the effect 
is particularly striking for small noise strengths F, and 
vanishes for F > 10"'^. 
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FIG. 3: Deflection angle -0 (see Fig. [ij and velocity v := |u| 
versus "dc-directionality" cj) (see fl}) for the same system as 
in Fig. Elbut with a = 7.4, and F = 4.4 ■ 10"* (solid), F = 
6.7-10"* (dashed). 



To better understand these findings we first focus on 
the deterministic dynamics (F = 0) in the simplest case 
when the ac-drive and the dc-bias are parallel and acting 
along one of the main symmetry axes of the periodic 
potential. Regarding the (1, 0) direction, i.e. (j) = a — 
0°, the lines y = n/2 constitute invariant sets of the 
deterministic dynamics (H]), stable for odd and unstable 
for even n. Thus, the particle motion is confined between 
two neighboring such lines and gencrically is attracted by 
the one with odd n for large times. Further, the velocity 
([2]) necessarily must follow the direction of the dc-force, 
i.e. "0 = 0°, and, in particular, vanishes for F — 0. This 
qualitative behavior remains unchanged in the presence 
of noise (F > 0). Analogous conclusions hold for the 
(0, 1) lattice direction. 

Turning to the (1,1) direction, i.e. = a = 45°, 
the lines x ^ y = n now constitute invariant sets of the 
deterministic dynamics ((T]) due to its invariance under 

51 : {x,y) I— > {y,x). Again, the particle motion must re- 
main confined between two adjacent such lines, implying 
for the velocity ([2]) that -0 = 0° and hence Vx = Vy. But 
now, the motion on the invariant lines may change its 
stability properties upon variation of a system parame- 
ter, and additional non-trivial attractors, not contained 
in any of the invariant lines, may arise. A typical exam- 
ple is shown in Fig. 21 We see that - depending on the 
driving amplitude a and possibly also on the seed r(0) 
- the orbit f{t) approaches a periodic or a chaotic long 
time behavior. The concomitant velocity ^ is still well 
defined but - in contrast to the noisy case F > - now 
may depend on the initial condition r(0). The "central" 
straight line in the bifurcation diagram for a < 3.22 be- 
longs to all the periodic attractors on the invariant sets 
X — y = n, exemplified by (a) in the upper inset of Fig. 
m and giving rise to a vanishing average velocity. At 
a 3.09 we observe the appearance of an additional pair 
of non-transporting [v = 0) periodic attractors, spon- 
taneously breaking 5*1 as well as the second symmetry 

52 ■ {x,y) I— *■ {—x,—y) of (HI), but still maintaining (up 
to translations) 5*1 o 5*2 : {x,y) 1— > {—y,—x), see (b) in 
Fig. m Thus, there are now three coexisting attractors 
within every unit cell of the periodic potential, one of 
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FIG. 4: Upper part: Bifurcation diagram for the unbiased 
{F = 0) deterministic (F = 0) dynamics ([T]) with a = 45°, 
SI = 3, and varying a. Shown is a stroboscopic representation 
of the attractors, governing the long-time behavior, by plot- 
ting x{jT) modulo L (L = 1), i.e. the reduced x-component 
at multiples j of the driving period T = 2n /Q,, for several dif- 
ferent seeds r(0) after initial transients have died out. Lower 
part: Corresponding a;-component of the average velocity ([2]) 
in units of L/T. Upper inset: Stable periodic orbits (attrac- 
tors) for a = 3.1 (a and b) and for a = 3.107 (c). Dotted: 
invariant sets x — y = n. Black "clouds" and "discs" repre- 
sent the potential U{r), corresponding to Fig. [T] after a 45° 
rotation. Lower inset: stable period-2 transporting orbit for 
a = 3.33. 



type (a), the others of type (b) and its image under 5*1. 
At a w 3.102 the pair of type (b) attractors exhibits 
a pitchfork bifurcation, which spontaneously breaks the 
Si o S2 symmetry (symmetry breaking bifurcation [3), 
resulting in four distinct, non-transporting attractors per 
unit cell. One of them is exemplified with (c) in Fig. 21 
its three "brothers" follow as mirror images with respect 
to the closest x ^ y — n and/or x + y = n lines. Upon 
further increasing a, a period doubling route to chaos fol- 
lows, which would be impossible without the preceding 
symmetry breaking bifurcation (l8| . The corresponding 
attractors loose stability by way of a crisis at a ~ 3.12. 
They reappear beyond a « 3.19 as "chaotic bands" , inter- 
rupted by "periodic windows" . Some of these windows 
exhibit attractors corresponding to phase-locked trans- 
porting orbits. The symmetry breaking bifurcation at 
a ~ 3.102 is pivotal for such transporting orbits: Since 
there is no systematic force {F = 0), which could fa- 



vor motion in one or the other direction, spontaneous 
symmetry breaking of all symmetries involving ^2 is an 
indispensable prerequisite for transporting orbits. The 
simplest and most prominent example arises within the 
periodic window at a « 3.33, exhibiting two attractors 
between any pair of adjacent x — y ~ n lines. One 
such orbit is exemplified in the lower inset of Fig. 21 
its "twin brother" follows as mirror image with respect 
to any x + y ~ n line. Apparently, this orbit arises by 
continuing the deformation of orbit (b), which leads to 
(c) even further, and rewiring one of its "arms" into the 
neighboring unit cell. The latter operation cannot be re- 
alized by a continuous deformation and hence one might 
guess that this somehow happens within the "gap" in the 
bifurcation diagram between a ~ 3.12 and a w 3.19. 

Due to the o S2 symmetry at = 0, oppositely 
transporting orbits co-exist and are stable within exactly 
the same range of the other system parameters. Applying 
a force F ^ 0, however, breaks the symmetry and hence 
the existence regions of the two orbits in parameter space 
no longer coincide. Closer inspection of how these regions 
change shape and size reveals that there are, for not too 
large \F\, parameter values, where the only stable orbit 
is the one which transports against the dc-bias F. In 
other words, we recover yet another example of "pure" 
ANM {tp = 180°) 0,[ll|, which furthermore turns out to 
survive even in the presence of (sufRciently weak) noise 
(F > 0), 

We finally address the case of arbitrary (but fixed) ori- 
entations of the external forcings. Without loss of gen- 
erality we restrict ourselves to < a < 90° but admit 
arbitrary cf). For such general driving directions the de- 
terministic dynamics ([1]) is not restricted any more by 
simple invariant sets. Rather, by means of extensive 
computer simulations we have found that transporting 
particle motion is created in a way very similar to the 
case exemplified above with Fig, 01 and typically "locks" 
to one of the three main symmetry axes (1,0), (0, 1) or 
(1, 1), depending on amplitude and frequency of the ac- 
drive. In particular, for F = 0, symmetry dictates the co- 
existence of transport into opposite directions for either 
of these basic orientations. Applying a (not too large) 
dc-bias Fe^ along a direction that is generally different 
from the one of the ac-drive, cf) ^ a, has two main effects: 
First, this co-existence is lifted, yielding parameter re- 
gions where only one transporting direction out of the 6 
different possible directions is stable. Second, new trans- 
port directions around the orientation of the bias force 
become accessible that follow a similar "locking scheme" 
as in [1^ upon variation of a system parameter. 

In other words, systematically changing, say, the driv- 
ing amplitude a leads to "jumps" in the deterministic 
transport direction where - due to the time-dependent 
ac-drive - also transport with velocity components oppo- 
site to the bias force occurs, in marked contrast to pre- 
vious findings in [l3|. The main effect of (weak) thermal 
noise is to "interpolate" between neighboring determin- 
istic directions, resulting in the smooth behavior shown 
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in Fig. [2] transport in virtually any direction but with 
emphasis on directions around the orientation of the dc- 
bias. Similarly, the variations of the deflection angle -ip 
observed in Fig. [3] result from a basically unchanged 
orientation of v into the "negative" (1,0) direction as 
long as the dc-bias has a non- vanishing component in the 
"positive" (1,0) direction (i.e. -90° < (j) < 90°) and a 
quick transition into the opposite situation when the dc- 
orientation (p moves into the complementary regime be- 
tween 90° and 270°. In any case, too large noise strengths 
r basically override the effects of the periodic potential 
and the system tends to return to the trivial behavior in 
the absence of the periodic potential. 

In conclusion, several quite astonishing linear response 
transport phenomena of a very simple and general non- 
equilibrium system have been observed and understood 



to the extent that an efficient and systematic search of 
pertinent parameter regions becomes easily feasible. The 
system is minimal in so far as any further reduction or 
simplification unavoidably rules out one of the indispens- 
able prerequisites, most notably the occurrence of spon- 
taneous symmetry breaking and chaos in the determin- 
istic limit. On the other hand, the effects are robust 
against a large variety of modifications and amendments 
of the model and hence should be realizable in several dif- 
ferent experimental systems [H, 0, 0, H, S [l^l for instance 
for particle sorting purposes. 
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